Inclusive dielectron production in proton-proton collisions at 2.2 GeV beam energy 
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Data on inclusive dielectron production are presented for the reaction p + p at 2.2 GeV measured 
with the High Acceptance DiElectron Spectrometer (HADES). Our results supplement data obtained 
earlier in this bombarding energy regime by DLS and HADES. The comparison with the 2.09 GeV 
DLS data is discussed. The reconstructed e^e~ distributions are confronted with simulated pair 
cocktails, revealing an excess yield at invariant masses around 0.5 GeV/c^. Inclusive cross sections 
of neutral pion and eta production are obtained. 



PACS numbers: 25.40Ep, 13.40.Hq 

I. INTRODUCTION 

The spectroscopy of e+e~ pairs offers a new approach 
to the study of baryon resonances excited in nucleon- 
nucleon reactions. Dilepton (that is e+e" or ob- 
servables provide indeed information on the electromag- 
netic structure of the resonances and, in the context of 
vector meson dominance, their coupling to the light vec- 



tor mesons Furthermore, dilepton spectroscopy al- 
lows to study the properties of hadrons produced and 
decayed in a strongly interacting medium. This is be- 
cause leptons do not themselves interact strongly when 
propagating through hadronic matter, that is, their kine- 
matics remains basically undistorted. For that reason 
they are used to probe medium modifications of hadrons 
intensively searched for in photon and proton-induced 
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reactions on nuclei as well as in heavy-ion collisions 
Transport models are commonly employed to de- 
scribe particle production and propagation through the 
medium, in particular when dealing with the complex 
dynamics of nucleus-nucleus reactions |3t-Q. The proper 
modeling of lepton pair production mechanisms requires 
a solid understanding of the underlying elementary pro- 
cesses, be it in terms of resonance excitations or in terms 
of a string fragmentation picture @ . 

The HADES experiment pursues a comprehensivepro- 
gram of dielectron emission studies in A*" -I- iV 0, IJ j in 
p + A , as well as in A -f A collisions (lol - [l^ . Inclusive 
e'^e~ production in p+p and p+d reactions had formerly 
been studied in the range of 1 - 5 GeV by the DLS exper- 
iment at the Bevalac [13], and more recently by HADES 
at 1.25 GeV [71] and 3.5 GeV 8]. In particular, the com- 
parison of the latter data sets with various model cal- 
culations demonstrated a need for improved theoretical 
descriptions. In this paper we supplement the available 
body of experimental results with data obtained on in- 
clusive e+e" production in the p+p — > p+p+e^ + e~ +X 
reaction at 2.2 GeV. A direct comparison with DLS data 
measured at 2.09 GeV [13] is presented. Furthermore, 
through the comparison with a calculated e"'"e~ cock- 
tail, we extract the inclusive production cross sections 
of tt" and rj mesons at 2.2 GeV. Our paper is organized 
as follows: Section 2 describes the experiment and the 
data analysis. In Sec. 3 e+e" pair spectra are presented 
and confronted with results from DLS. In Sec. 4 the pair 
spectra are compared to calculated dielectron cocktails. 
In Sec. 5 we discuss inclusive meson production cross sec- 
tions and, finally, in Sec. 6 we summarize our findings. 



II. THE EXPERIMENT 

The six-sector High- Acceptance DiElectron Spectrom- 
eter (HADES) operates at the GSI Helmholtzzentrum fiir 
Schwerionenforschung in Darmstadt taking beams from 
the heavy-ion synchrotron SIS18. Technical aspects of 
the detector are described in [T^. Its main component 
serving for electron and positron selection is a hadron- 
blind Ring-Imaging Cherenkov detector (RICH). Further 
particle identification power is provided by the time of 
flight measured in a plastic scintillator wall (TOE), the 
electromagnetic shower characteristics observed in a pre- 
shower detector, and the energy- loss signals from the 
scintillators of the TOE wall. 

In the experiment discussed here [H] a proton beam 
with a kinetic energy of Tp = 2.2 GeV (corresponding 
to a cm. energy ^^snn = 2.765 GeV) and an intensity 
of about 10^ particles per second impinged on a 5 cm 
long liquid hydrogen cell with a total areal thickness of 
0.35 g/cm^. The online event selection was done in two 
steps: (1) a l**-level trigger (LVLl) selected events with 
an overall multiplicity of at least four charged hits in the 
TOE wall with additional topological conditions (two op- 
posite sectors hit, two hits at polar angles < 45°), and 



(2) a 2"''-level trigger (LVL2) required an electron or 
positron candidate. This trigger scheme was in fact op- 
timized for measuring exclusive e+e^ production in the 
p + p -J_p + P + V reaction with a subsequent 77 Dalitz 
decay [16[. Note that such a trigger condition still al- 
lows to study, albeit with a bias, inclusive e"^e^ emis- 
sion. Indeed, because of overall charge conservation in 
the p + p reaction, the dielectron is always accompanied 
by at least two more charged particles. The resulting 
trigger bias has been studied in simulations as a function 
of various pair observables, in particular pair mass and 
pair transverse momentum, providing a correction as well 
as an estimate of the resulting systematic error (of order 
20%). Eor normalization purposes, p + p elastic scatter- 
ing events were recorded concurrently with an additional 
scaled-down (by a factor 32) LVLl trigger condition re- 
quiring only two charged hits in opposite HADES sec- 
tors. Thus, in total 2.7 x 10^ LVLl events were recorded, 
4.1 X 10^ of which fulfilled the LVL2 condition. 

Dielectrons (that is e+e^ pairs) were reconstructed 
following the procedures described in detail in [ll|, [ij] : 
(1) leptons were identified based on various detector ob- 
servables, (2) an efficiency correction was applied, (3) 
opposite-sign leptons were combined into pairs, (4) the 
background of uncorrelated (and partially correlated) 
pairs representing the combinatorial background (CB) 
was subtracted using the same-event geometric mean of 
like-sign pairs, (5) the correction for the LVLl trigger 
bias was applied, and finally (6) the resulting inclusive 
e^e~ distributions were normalized to the reconstructed 
yield of elastically scattered protons into the HADES ge- 
ometric acceptance (see [1^ for details) . As no dedicated 
start detector was present in this experimental run, the 
start time for the time-of-fiight measurement was recon- 
structed event-by-event from the most optimal fit of dif- 
ferent event hypotheses to the global event data p^ . 



III. RESULTS 

A. Invariant mass spectra 

Figure [T] shows the differential e~^e~ cross section 
da/dMee obtained after correcting the reconstructed pair 
yield for efficiency, combinatorial background, and trig- 
ger bias. As explained above, the absolute normalization 
was done using the known p + p elastic scattering cross 
section with the help of the concurrently measured yield 
of elastic events [16'] . The data are presented with anal- 
ysis cuts on single-lepton momentum, p^ > 0.1 GeV/c, 
and on pair opening angle, ^ee > 9°. The total num- 
ber of e+e^ signal pairs contributing to this spectrum is 
around 19,000 from which 2,000 pairs are located above 
the tt" Dalitz region {Mee > 0.15 GeV/c^). To illustrate 
the significance of the reconstructed dielectron signal, the 
signal-ovcr-CB ratio is also shown as an inset in Fig. [T] 
Note that the kinematic cutoff corresponding to the beam 
energy of 2.2 GeV is at a pair mass of 0.89 GeV/c^. 
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FIG. 1: (Color online) Differential e^e~ cross section 
do/dMee measured in the 2.2 GeV p + p reaction within the 
HADES acceptance (including pe > 0.1 GeV/c and 6ee > 9° 
cuts). The data are efficiency corrected and CB subtracted; 
the insert shows the signal-over-CB ratio (S/CB). Point-to- 
point statistical and systematic errors are indicated by (black) 
vertical bars and (red) horizontal ticks, respectively. 

The obtained cross section is combined in Fig. [5] with 
inclusive data from other HADES p + p runs done at 
1.25 [7] and 3.5 GeV Q bombarding energy, respectively. 
While the 3.5 GeV data were recorded in the same detec- 
tor acceptance as the present 2.2 GeV data, the 1.25 GeV 
data were adjusted for differences in the magnetic field 
strengths and analysis cuts used. This way the three data 
sets are compared within the same instrumental accep- 
tance, revealing the strong beam energy dependence of 
dielectron production, particularly at large pair masses. 
Note, however, that between the three beam energies 
the average momentum of the involved single lepton dis- 
tributions differs significantly, resulting in substantially 
different fractions of accepted pairs, particularly for low 
masses. 



B. Comparison with DLS 

As the former DLS experiment provided e^e^ data for 
the p + p reaction at 2.09 GeV we can make a direct 
comparison along the line already used to confront the 
DLS and HADES C-l-C data obtained at I A GeV 10, 17]. 
As the HADES geometrical acceptance is substantially 
broader than the DLS one such a comparison can be done 
by projecting the reconstructed HADES dielectron yields 
d^N/dMeedP±dY through the DLS acceptance filter [13] 
(here P± and Y are the e^e~ pair transverse momentum 



FIG. 2: (Color online) Systematics of e'^e differential pro- 
duction cross sections da/dMes measured in p + p reactions 
at 1.25 GeV (squares), 2.2 GeV (circles), and 3.5 GeV (trian- 
gles), all obtained within the HADES acceptance, efficiency 
corrected, and CB subtracted (including pe > 0.1 GeV/c and 
Oee > 9° cuts). The 1.25 GeV data, taken from 0], are ad- 
justed for the present, more restrictive detector acceptance 
(i.e. stronger magnetic field and explicit lepton momentum 
cut of 0.1 GeV/c); the 3.5 GeV data are taken from 8]. Only 
statistical error bars are shown. 

and rapidity, respectively). In fact, as DLS applied to 
their p + p data additional cleaning cuts [H, [3] - 0.1 < 
Mee < 1-25 GeV/c2, Pl < 1.2 GeV/c, 0.5 < V < 1.7, 
and 6*6 > 21.5° - an extrapolation of the HADES yield to 
rapidities above 1.9, as applied in [lOj, is not needed here. 
The result of this filtering procedure is shown in Fig.[3Ua) 
for the pair mass distributions da/dM^e and in (b) for the 
pair transverse momentum spectrum l/(27rPx) da/dP±, 
the latter one with the condition Mee > 0.15 GcV/c^. 

It is apparent that, within statistical and systematic 
uncertainties, the HADES and the DLS data are in good 
agreement. This result suggests that our result, together 
with the data obtained by the DLS energy scan [l3], 
can be used to constrain the various models aimed at 
describing dielectron production in few-GeV elementary 
reactions [8j]- 



IV. COMPARISON WITH A SIMULATED 
COCKTAIL 

The experimental pair distributions are next compared 
to a calculated e+e^ cocktail. For this, p + p reactions 
were simulated with the event generator Pluto [l^, [2^ 
and filtered through the HADES acceptance. The sim- 
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FIG. 3: (Color online) Direct comparison within the DLS 
acceptance (see text for details) of the e^e~ cross sections 
measured by HADES in p + p at 2.2 GeV (closed circles) and 
by DLS at 2.09 GeV (open crosses, taken from [l3|). The 
pair mass distributions (a) and pair transverse-momentum 
distributions (b) are confronted. Error bars are statistical 
only; additional systematic errors (not shown) are 23% for 
DLS jj] and 29% for HADES. 



Illation included the following e+e pair sources: (1) 



TfO ^ 7e+e-, (2) ^ 76+6", (3) A(1232) -> A^e+e", 
(4) uj TT^e+e", (5) w e+e", and (6) p° e+e~ with 
dielectron branching ratios of mesons taken from [2l| and 
that of the A(1232) as calculated in At the present 
bombarding energy the production of tt" and 77 mesons 
is known to proceed mostly via resonance excitation 
(e.g. R = A(1232), Ar*(1440),iV*(1520), A^*(1535), etc.) 
and is dominated by one-meson and two-meson channels 



For and uj production we have, however, 
assumed phase-space population in the pp — >■ ppX reac- 
tions {X — or w) with no further attempt at a more 
refined description of the high-mass region. Note that 
some of the excited resonances i?, mostly the A, con- 
tribute also directly to the dielectron yield via their elec- 
tromagnetic Dalitz decay R Ne~^e~. In our cocktail 
calculation we have only taken into account the A'^ and 
A+ contributions following the prescription from Ref. Q . 
The production cross sections used in the simulation were 
taken as follows: 

1. Inclusive 7r° production ~ 14 mb (adding all ob- 
served inelastic channels contributing to tt*^ pro- 
duction from Ref. [25| gives a lower limit of about 
12 mb, whereas 14 mb are needed to fully exhaust 
the measured Dalitz yield). 

2. Inclusive i] production - in the range of 0.26 - 
0.35 mb; a lower limit of 0.14 mb is given by the 
known exclusive 77 production [TgI . [26j . 

3. Vector meson production - 0.01 mb exclusive uj pro- 
duction [l^] and assuming likewise for the but 
no (f) contribution. 

4. Inclusive A'^'+(1232) excitation - in the range 10 - 
21 mb. 

The extremes of the cross section range used for A pro- 
duction correspond to the following two scenarios: (I) 
assumes that pion production is mediated completely by 
single A excitation only, that is cta — 3/2cr^o (from 
isospin addition rules), resulting in the upper value of 
21 mb; scenario (II) sums explicitly the A contributions 
from one-pion and two-pion production channels. In the 
latter the one-pion part of 3.6 mb is taken from a res- 
onance model fit to exclusive pion production data [l6l | 
whereas the two-pion part of 6.4 mb is taken from the 
effective Lagrangian model of two-pion production pre- 
sented in Ref. j24| . Under the assumption and as sug- 
gested indeed by various calculations that dielec- 
tron production is dominated in the mass range 0.15 - 
0.45 GeV/c^ by the A and contributions, the total 
yield measured for these masses can be used to constrain 
the T] contribution for any assumed A cross section. In 
other words, the 77 and A contributions are complemen- 
tary. The extracted rj cross section will evidently have a 
model dependence which, however, turns out to be small. 

The resulting e+e^ cocktails, filtered with the HADES 
acceptance, are overlayed in Figs. |4] and [5] with the data. 
Up to masses of ~ 0.45 GeV/c^ the agreement is good in 
both observables, M^e and Pi, although at higher masses 
the measured yield is not matched. Integrating up to 
0.15 GeV/c^ the low-mass region, dominated evidently 
by the 7r° Dalitz contribution, and correcting for the de- 
tector acceptance, allows to fix the inclusive 7r° produc- 
tion cross section at cr^o = 14 ±3.5 mb. The quoted 25% 
error is determined mostly by systematic effects (normal- 
ization, trigger bias correction, acceptance correction). 
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In a similar way, the integrated yield from the mass re- 
gion of 0.15 - 0.45 GeV/c^ has been used to extract the 
cross section of inclusive rj production after correcting, 
as stated above, for the A contribution. The range of 
assumed A cross sections used in the simulation leads 
consequently to a corresponding range of rj production 
cross sections, indicated by the hatched bands shown in 
the two figures. In scenario (I), where all tt'^ production 
goes through A excitation and decay, cta = 21 mb and 
a„ = 0.26 mb. In scenario (II), based on the model of 
[2J|, various nucleon resonances contribute to pion pro- 
duction, such that CTA = 10 mb and cr^ = 0.35 mb. We 
feel that this model dependence is relatively small and 
propose to use the average of the two results, namely 
cr^ — 0.31 ± 0.08 ± 0.05 mb, where the first error is again 
mostly ruled by systematic effects (normalization, trig- 
ger bias and acceptance corrections) while the second one 
covers the model dependence. 

In the pair mass range 0.45 - 0.60 GeV/c^ our simu- 
lation underestimates grossly the observed yield. This 
is also visible in the comparison of transverse momen- 
tum distributions shown on Fig. [5] Clearly additional 
dielectron sources are needed, among which one has to 
consider the decays of N* resonances, e.g. the Af*(1520) 
and A^*(1720) known to couple strongly to the p, as well 
as a possible general enhancement due to vector meson 
dominance form factors of the nucleon resonances y . 




V. INCLUSIVE MESON PRODUCTION 

The inclusive tt*^ and rj production cross sections ob- 
tained in the present analysis can be combined with the 
result from our p + p runs at 1.25 and 3.5 GeV to in- 
vestigate the excitation function of meson production in 
the few-GeV regime. Figure H] shows these cross sections 
as function of ^^snn together with exclusive data. A 
wealth of information on exclusive pion production in 
nucleon-nucleon reactions has indeed been accumulated 
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over the past 50 years (see [22, l23| for reviews and [25 1 
for a compilation). Fits to exclusive tt" production cross 
sections in one-pion and two-pion channels from [22j | are 
shown in Fig. [S] These processes are quite well under- 
stood in terms of nucleon resonance excitations within 
the resonance models [13, 113] . Data on tt" production in- 
volving three or more pions in the final state is, however, 
still scarce and incomplete, although such processes can 
be expected to contribute substantially at beam energies 
above a few GeV. Indeed, from an extrapolation of three- 
pion production data to Tp = 2.2 GeV published recently 
(28{ it is estimated that the contributions of the tt+tt^tt'^ 
and TT^TT^TT^ channels could add up to as much as 0.5 mb. 
Figure |6] shows in fact very clearly that around 2.2 GeV 
bombarding energy inclusive pion production stops to be 
fully exhausted by the sum of one- and two-pion channels 
only. 

Turning to rj production we can do a similar compari- 
son. Here, data has again been only available for the ex- 
clusive channel. A compilation of exclusive rj production 
cross sections (from [2^, extended with recent HADES 
results n^) is depicted in Fig. |6l as well as the corre- 
sponding resonance model calculation of Teis et al. [29{ . 
Assuming that t] production is mediated solely by the 
iV* (1535) resonance, the latter gives a reasonable descrip- 
tion of the data. Just like in the case of pion production, 
our inclusive cross sections largely exceed the exclusive 
ones in the energy range investigated here. Because a mi- 
croscopic description of multi-particle production is not 
yet at hand, transport models often make use of cross 
section parameterizations based on the Lund string frag- 
mentation model (LSM) Sibirtsev's parameteriza- 
tion of rj production fSO'l , based on the LSM and shown 
as dot-dashed line in Fig. |6l turns out to be in reason- 
able agreement with our inclusive result, although the 
intended validity range of the LSM is in fact at much 
higher beam energies. 



FIG. 4: (Color online) Pair mass distribution da/dM^e mea- 
sured in 2.2 GeV p + p reactions (full circles) compared with 
simulated Pluto [19l.[2(]|| cocktails of dielectron sources filtered 
through the HADES acceptance. The shaded bands delimit 
the range of modeled A and r) contributions - as discussed in 
the text - with the dashed delimiters corresponding to sce- 
nario (I) and the solid ones to scenario (II). Only statistical 
errors are shown. 



VI. SUMMARY 

To summarize, we have presented data on inclusive 
e'''e^ production in the reaction p + p td, 2.2 GeV beam 
kinetic energy. The measured dielectron cross sections 
are in good agreement with the DLS result obtained ear- 
lier at 2.09 GeV. The employed cocktail of e+e~ sources 
does not saturate our data at invariant mass around 
0.55 GeV/c^. Furthermore, inclusive tt" and rj produc- 
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FIG. 5: (Color online) Pair transverse momentum distributions da/dP± measured in 2.2 GeV p+p reactions within the HADES 
acceptance. Three mass bins are shown: (a) Mee < 0.15, (b) 0.15 < Ale^ < 0.45, and (c) Mee > 0.45 GeV/c^. The simulated 
Pluto cocktails are shown as well, with line styles as in Fig. [J] 
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FIG. 6: (Color online) tv" and rj production cross sections 
in p + p reactions as a function of the cm. energy yijvjv 
(bottom scale) and beam kinetic energy Tp (upper scale). The 
present inclusive results are shown as full triangles and circles, 
respectively, together with more HADES data obtained at 
bombarding energies of 1.25 GeV 0] and 3.5 GeV Fits 
to a compilation of Itv and Ivr + 27r cross sections [23 ] are 
shown as solid and long-dashed curves, respectively. Open 
squares are ri exclusive cross sections taken from j2i_j23|; the 
solid square is the exclusive HADES point from [la ]. The 
short-dashed curve corresponds to the resonance model [29l ]. 
and the dot-dashed curve is the parametrization of inclusive 
■q production from [30l ]. 



tion coss sections were deduced, extending the world 
body of meson production data. Such data, besides hav- 
ing their own interest in the context of medium-energy 
hadron reactions, represent valuable input for the analy- 
sis and simulation of proton-nucleus and heavy-ion colli- 
sions in the same energy regime. 

Acknowledgments 



The collaboration gratefully acknowledges the sup- 
port by CNRS/IN2P3 and IPN Orsay (France), by SIP 
JUG Cracow (Poland) (NN202 286038, NN202198639), 
by HZDR, Dresden (Germany) (BMBF 06DR9059D), 
by TU Miinchen, Garching (Germany) (MLL Miinchen, 
DFG EGlust 153, VH-NG-330, BMBF 06MT9156 TP5, 
GSI TMKrue 1012), by Goethe-University, Frankfurt 
(Germany) (HA216/EMMI, HIC for FAIR (LOEWE), 
BMBF 06FY9100I, GSI F&E), by INFN (Italy), by 
NPI AS CR, Rez (Czech Repubhc) (MSMT LG07050, 
GAASCR IAA100480803), by USC - Santiago de Gom- 
postela (Spain) (CPAN:GSD2007-00042). 



[1] M I. Krivoruchenko, B V. Martemyanov, A. Faessler, and 

C. Fuchs, Ann. Phys. 296, 299 (2002). 
[2] S. Leupold, V. Metag, and U. Mosel, Int. J. Mod. Phys. 

E 19, 147 (2010). 
[3] E. L. Bratkovskaya and W. Gassing, Nucl. Phys. A 807, 

214 (2008). 



J4] K. Schmidt et al., Phys. R ev. C 79, 064908 (2009). 
J5] O. Buss et al., arXiv: 1106. 134^1 [hep-ph], submitted to 
Phys. Rep. 

J6] B. Nilsson-Almqvist and E. Stenlund, Comp. Phys. 

Comm. 43, 387 (1987). 
J7] G. Agakishiev et al. (HADES Collaboration), Phys. Lett. 



7 



B 690, 118 (2010). 
[8] G. Agakishiev et al. (HADES Collaboration), 
larXiv:1112.3607t ^2 [nucl-ex], submitted to Eur. Phys. J. 
A. 

[9] M. Weber et al. (HADES Collaboration), J. Phys. Conf. 

Ser. 316, 012007 (2011). 
[10] G. Agakishiev et al. (HADES Collaboration), Phys. Lett. 

B 663, 43 (2008). 
[11] G. Agakishiev et al. (HADES Collaboration), Phys. Rev. 

Lett. 98, 052302 (2007). 
[12] G. Agakishiev et al. (HADES Collaboration), Phys. Rev. 

C 84, 014902 (2011). 
[13] W. Wilson et al. (DLS Collaboration), Phys. Rev. C 57, 

1865 (1998). 

[14] G. Agakishiev et al. (HADES Collaboration), Eur. Phys. 

J. A 41, 243 (2009). 
[15] B. Sailer, doctoral thesis. Technical University, Miinchen 

(2007). 

[16] G. Agakishiev et al. (HADES Collaboration), 
larXiv:1203.1333t rl [nucl-ex], submitted to Eur. Phys. J. 
A. 

[17] R. J. Porter et al. (DLS Collaboration), Phys. Rev. Lett. 

79, 1229 (1997). 
[18] H. Matis, LBNL, private communication. 
[19] I. Frohlich et al., Proceedings of the XI Interna- 



tional Workshop on Advanced Computing and Analy- 
sis Techniques, Amsterdam (The Netherlands) 2007, PoS 
(ACAT) 076. 

[20] F. Dohrmann et al., Eur. Phys. J. A 45, 401 (2010). 
[21] K. Nakamura et al. (Particle Data Group), J. Phys. G 

37, 075021 (2010). 
[22] J. Bystricky et al., J. Physique 48, 1901 (1987). 
[23] P. Moskal, M. Wolke, A. Khoukaz, and W. Oelert, Prog. 

Part. Nucl. Phys. 49, 1 (2002) and refs. therein. 
[24] X. Cao, B.-S. Zou, and H.-S. Xu, Phys. Rev. C 81, 065201 

(2010). 

[25] A Baldini, V. Flaminio, W. G. Moorhead, and 

D. R. O. Morrison, in Landolt-Bornstein, New Series 

I/12b, (Springer 1988). 
[26] F. Balestra et al. (DISTO Colaboration) , Phys. Rev. C 

69, 064003 (2004). 
[27] M. Abdel-Bary et al. (COSY-TOF Collaboration), Phys. 

Lett. B 647, 351 (2007). 
[28] C. Pauly et al. (CELSIUS- WASA Collaboration), Phys. 

Lett. B 649, 122 (2007). 
[29] S. Teis et al., Z. Phys. A 356, 421 (1997). 
[30] A. Sibirtsev, W. Gassing, and U. Mosel, Z. Phys. A 358, 

357 (1997). 



